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METHODS FOR DETERMINING THE TRUE 
SIGNAL OF AN ANALYTE 

This application is based on, and claims the 
benefit of, U.S. Provisional Application No. 60/248,259, 
5 filed November 14, 2000, entitled Testing for 

Differentially-Expressed Genes by Maximum Likelihood 
Analysis of Microarray Data and claims benefit of, U.S. 
Provisional Application No. 60/266,388, filed February 2, 
2001, entitled Methods for Determining the True Signal of 
10 an Analyte, which are incorporated herein by reference. 

This invention was made with government support 
under grant number T32 HG 000-35 awarded by the National 
Institutes of Health and grant number 
DE-FG03-98ER62652/A000 awarded by the United States 
15 Department of Energy. The United States Government has 
certain rights in this invention. 

BACKGROUND OF THE INVENTION 

The invention relates generally to quantitative 
expression analysis, and more particularly, to methods 
20 for identifying significant differences in gene 
expression. 

Although all cells in the human body contain 
the same genetic material, the same genes are not active 
in all of those cells. Alterations in gene expression 
25 patterns or in a DNA sequence can have profound effects 
on biological functions. These variations in gene 
expression are at the core of altered physiologic and 
pathologic processes. In the past, determinations of 
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differential gene expression only focused on a few genes 
at a time. DNA microarrays, devices that consist of 
thousands of immobilized DNA sequences present on a 
miniaturized surface, have revolutionized the study of 
5 gene expression and are now a staple of biological 
inquiry into gene expression and genetic variations. 
Arrays are used to analyze a sample for genotyping or for 
patterns of gene expression. Using the microarray, it is 
possible to observe the expression level changes in tens 
10 of thousands of genes over multiple conditions, all in a 
single experiment. Depending on the conditions assayed, 
differentially-expressed genes may be implicated in 
cancer, aging, or a metabolic pathway of interest. 

Generally, microarrays are prepared by binding 

15 DNA sequences to a surface such as a nylon membrane or 
glass slide at precisely defined locations on a grid. 
Using an alternate method, some arrays are produced using 
laser lithographic processes and are referred to as 
biochips or gene chips. For genotyping analysis, the 

20 sample is genomic DNA. For expression analysis, the 

sample is cDNA, DNA copies of mRNA. The DNA samples are 
tagged with a radioactive or fluorescent label and 
applied to the array. Single stranded DNA will bind to a 
complementary strand of DNA. At positions on the array 

25 where the immobilized DNA recognizes a complementary DNA 
in the sample, binding or hybridization occurs. The 
labeled sample DNA marks the exact positions on the array 
where binding occurs, allowing automatic detection. The 
output consists of a list of hybridization events, 

30 indicating the presence or the relative abundance of 
specific DNA sequences that are present in the sample. 



DNA array technology provides a method for rapid 
genotyping, facilitating the diagnosis of diseases for 
which a gene mutation has been identified as well as for 
diseases for which known gene expression biomarkers of a 
5 pathologic state, or signature genes, exist. 

A crucial step in the analysis of expression 
data is determining which genes are expressed differently 
between two cell populations. Usually, a gene is said to 
be "differentially-expressed" if its ratio of expression 

10 level in one population to expression level in a second 
population exceeds a certain threshold. This threshold 
is set based on the observation that in control 
experiments where the two cell populations are identical, 
few if any genes have expression ratios exceeding the 

15 threshold. However, it is common knowledge that this 
approach is imprecise, because the uncertainty in the 
expression ratio is greater for genes that are expressed 
at low levels than for those that are highly expressed. 
More sensitive methods have been employed in a few cases, 

20 but development of a general, formal statistical test for 
identifying differentially-expressed genes has remained 
an open problem. 

Thus, there exists a need for a mathematical 
model of the variability observed over repeated 
25 observations of intensities for biomolecules represented 
on an array. The present invention satisfies this need 
and provides related advantages as well. 
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SUMMARY OF THE INVENTION 

The invention relates to a method of 
determining a true signal of an analyte, comprising (a) 
measuring an observed signal x for one or more analytes, 
5 and (b) determining a mean signal (p.) and a system 
parameter (p) for said analyte that produce enhanced 
values for a probability likelihood of said observed 
signal, said observed signal being related to said mean 
signal by an additive error (5) and a multiplicative 
10 error (£), wherein said system parameter specifies 
properties of said additive error (5) and said 
multiplicative error (e). 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows the (A) increase of standard 
15 deviation; (B) increase of correlation with absolute 

level of intensity x' or y*; and (C) normal probability 
plot for the 80 samples of x 1 pertaining to a single, 
representative gene. 

Figure 2 shows scatter plots of estimated u y 
20 versus u x for each gene represented on the whole-yeast 
genome microarray, for (A) the control experiment YPRG 
versus YPRG and (B) the YPR versus YPRG comparison, while 
(C) shows the distribution of four (x,y) pairs for two 
genes in the YPR versus YPRG comparison . 

25 Figure 3 shows array images corresponding to 

hybridizations performed for each of eight controlled 
GAL80 ratios where four (x,y) intensity measurements per 
gene were obtained at each controlled GAL80 ratio by 
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using (A) two spots from a forward Cy3:Cy5 labeling 
scheme and (B) two spots from a reverse Cy5:Cy3 labeling 
scheme and (C) comparison of each controlled ratio to 
measured ratio (y/x) for the forward-array (red dots) or 
5 reverse-array (green dots) . 

DETAILED DESCRIPTION OF THE INVENTION 

The invention provides a method of determining 
relative amounts of an analyte between samples. The 

10 invention also provides a method of determining the true 
signal of an analyte. The method of the invention 
accounts for multiplicative and additive errors 
influencing the observed signals for an analyte and 
estimates system parameters based on the observed signals 

15 using maximum likelihood estimation. By presenting an 

error model and associated significance test, the methods 
of the invention provide a substantial improvement over 
current thresholding schemes. One advantage of the error 
model is that the system parameters inherently specify 

20 the properties of both the additive and multiplicative 
error terms. The method of the invention further 
provides for the performance of a generalized likelihood 
ratio test for each analyte to determine whether the 
amounts are relatively different. 

25 In one embodiment, the method of the invention 

provides a refined test for comparison of differentially 
expressed genes that does not rely on gene expression 
ratios, but directly compares a series of repeated 
measurements of two observed intensities for each gene. 

30 In this regard, the method of the invention utilizes an 
error model and an associated significance test to 
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determine whether the observed amounts of genes are 
significantly different between the two or more 
conditions being compared. 

As used herein, the term "analyte" refers to a 
5 molecule whose presence is measured. An analyte molecule 
can be essentially any molecule for which a detectable 
probe or assay exists or can be produced by one skilled 
in the art. For example, an analyte can be a 
macromolecule such as a nucleic acid, polypeptide or 

10 carbohydrate, or a small organic compound. Measurement 
can be quantitative or qualitative. An analyte can be 
part of a sample that contains other components or can be 
the sole or major component of the sample. Therefore, an 
analyte can be a component of a whole cell or tissue, a 

15 cell or tissue extract, a fractionated lysate thereof or 
a substantially purified molecule. Moreover, an analyte 
can incorporate a second molecule, for example, a 
detectable moiety such as a dye, radiolabel, heavy atom 
label, or other mass label, a f luorochrome, a 

20 ferromagnetic substance, a luminescent tag or a 

detectable binding agent such as biotin. The analyte can 
be attached in solution or solid-phase, including, for 
example, to a solid surface such as a chip, microarray or 
bead. 

25 As used herein, the term M sample" refers to the 

substance containing the analyte. It can be 
heterogeneous or homogeneous. Examples of heterogeneous 
samples include tissues, cells, lysates and fractionated 
portions thereof. Homogeneous samples include, for 

30 example, isolated populations of polypeptides, nucleic 
acids or carbohydrates. A sample can also be a purified 
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analyte, free from like or non-like molecules. All of 
such substances are included within the meaning of the 
term so long as the substance contains the analyte. In 
addition to containing the analyte, a sample further can 
5 contain one or more additional components such as a 
buffer, detectable moiety, nucleic acids, polypeptides, 
carbohydrates or any other substance or molecule. 

As used herein, the term "signal" is intended 
to mean a detectable, physical quantity or impulse by 

10 which information on the presence of an analyte can be 
determined. Therefore, a signal is the read-out or 
measurable component of detection. A signal includes, 
for example, fluorescence, luminescence, colorimetric, 
density, image, sound, voltage, current, magnetic field 

15 and mass. Therefore, the term "observed signal," as used 
herein is intended to mean the actual quantity detected 
of the measured analyte in a particular detection system. 
An observed signal can include subtraction of 
non-specific noise. An observed signal can also include, 

20 for example, treatment of the measured quantity by 

routine data analysis and statistical procedures which 
allow meaningful comparison and analysis of the observed 
values. Such procedures include, for example, 
normalization for direct comparison of values having 

25 different scales, and filtering for removal of aberrant 
or artif actual values. A "mean signal" as used herein, 
refers to the true or inherent quantity of the measured 
analyte. A mean signal therefore corresponds to the 
detectable quantity of the analyte independent of 

30 variation in the assay or detection system. 
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As used herein, the term "sample pairs" refers 
to two samples containing analytes to be compared. The 
analytes to be compared within the two samples can be 
different, or they can be substantially the same species 
5 of analyte but subjected to distinct conditions or 
obtained from distinct sources. Therefore, the term 
"mean signal pairs/' as used herein, refers to the two 
true signals, one per analyte, associated with a sample 
pair. Similarly, when more than two analytes are being 
10 compared, the terms "sample sets" and "mean signal sets" 
are intended by analogy to reference the multiple samples 
containing the analytes and the corresponding multiple 
true signals, respectively. 

As used herein, the term "system parameter" 
15 refers to the properties of the noise of the system, such 
as non-analyte, non-specific background signals. 
Therefore, the system parameter, designated p, is a 
measure of the error of the system and corresponds to 
undesirable or interfering signals that distort the true 
20 signal. 

As used herein, the term "significantly 
unequal" refers to two analytes that have a meaningful 
difference in signal. Therefore, significantly unequal 
signals refers to two or more signals whose difference is 
25 caused by something other than chance, including 
variation or error in the system. 

The invention provides a method of determining 
a true signal of an analyte. The method consists of 
measuring an observed signal x for one or more analytes 
30 and determining a mean signal (u) and a system parameter 
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(3) for the analyte that produce enhanced values for the 
probability likelihood of the observed signal, which is 
related to the mean signal by an additive error (5) and a 
multiplicative error (e), where the system parameter 
5 specifies the properties of the additive error (5) and of 
the multiplicative error ( 2 ) . 

The invention further provides a method 
determining relative amounts of an analyte between 
samples- The method consists of measuring observed 

10 signals x and y for an analyte within two or more sample 
pairs, determining a mean signal pair per analyte (u) and 
a system parameter ((3) for each sample pair, that produce 
enhanced values for the probability likelihood of the 
observed signals, which are related to the mean signal by 

15 an additive error (5) and a multiplicative error (£), 
where the system parameter specifies the properties of 
the additive error (5) and the multiplicative error (£). 

The methods of the invention permit 
determination of the mean signal, which is the true 

20 amount of an analyte, by taking into account both 

multiplicative and additive error contributions to each 
observed signal. The methods of the invention further 
allow accurate determination of relative amounts of an 
analyte between samples. A maximum-likelihood approach 

25 is used to fit the model to observed signals of the 
analyte. The method of the invention can be used to 
monitor error introduced by intrinsic or extrinsic 
factors, to monitor total amount of error over time as 
well as to isolate or identify particular samples that 

30 have a higher error than normally observed. Therefore, 
the methods of the invention can be used to detect error 
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introduced during any step in the analyte preparation and 
measurement. Additionally, the methods of the invention 
can be used, for example, to detect total error of the 
system or to separate and dissect biological or other 
5 intrinsic sample error from assay and procedure error* 
Thus, the methods of the invention allow quantitative 
analysis of the mean or true amount of an analyte at any 
given end point in a procedure as well as allow 
dissection of the system or procedure to quantitatively 
10 determine either or both intrinsic or extrinsic error 
introduced at any given step of the procedure. 

Likelihood methods use statistical data and 
probability models to provide optimal use of statistical 
information. Because likelihood methods provide a 
specific description of the pattern of variation in data, 
these methods can be used for estimation and hypothesis 
testing, which is a formal process of using data to make 
statistically meaningful decisions such as whether 
relative amounts of analyte are significantly different 
between samples. Therefore, the methods of the invention 
determine, by formal estimation procedures, the mean 
signal of an analyte or a comparison of mean signals to 
provide the relative levels of the corresponding analyte. 
The comparison of mean signals can be for the same 
analyte subjected to two or more different conditions, 
different analytes under the same conditions or any 
combination thereof. 

For comparison of two signals, the 
maximum-likelihood approach provided by the invention has 
30 several advantages over currently accepted ratio-based 
significance tests. In the ratio-based method, the 
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expression ratio for the two signals to be compared is 
computed and compared to a control or reference ratio. 
For example, where the relative level of an analyte is to 
be compared under two different conditions, the ratio 
5 r^Xi/yi is computed for analyte i for the two conditions 
x and y, and compared to a reference ratio of known 
analyte signals. A ratio that differs from the reference 
ratio, for example, as r t >r c or r i <l/r c identifies the 
analyte levels under the two conditions as being 

10 meaningfully different. This ratio-based method has been 
widely used in fields that compare, for example, the 
differences in expression of RNA or protein under two 
different conditions. The method has been particularly 
applicable to large scale expression analysis such as 

15 those utilizing microarray formats. However, the 

ratio-based method for statistical analysis of signal 
data combines observed signals into a single ratio, which 
necessarily results in the loss of absolute signal 
information. Moreover, when repeated samples per analyte 

20 are available, common practice is to compute the ratio of 
averaged signals, again discarding useful information. 

The methods of the invention are generally 
applicable to measure any analyte that serves as a sample 
or is contained in a sample so as to allow for detection 
25 of the presence of the analyte. As will be described in 
further detail below, detection of the analyte signal can 
be by any means as long the observed signal allows for 
determination of a mean. 

Once a signal indicating the presence of an 
30 analyte has been observed, the methods of the invention 
can be used to determine the true or mean signal of the 
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analyte. The true signal of an analyte is independent of 
experimental variation or error introduced prior to or 
during detection of the observed signal. Removal of such 
error in a signal allows for more accurate quantitation 
5 of an analyte and reproducibility of measurements. 
Therefore, the true or mean signal of an analyte is a 
measurement of the true or actual level of that analyte. 
Moreover, through determination of the true signal, the 
methods of the invention can measure the reproducibility 
10 of steps in a process such as, for example, manipulations 
prior to the determination of the observed signal. 

The methods of the invention are applicable to 
the measurement of analytes and determination of true 
signals in both biological and non-biological settings. 

15 For example, in a biological setting, experimental error 
can be classified into at least two categories. 
Biological error is one such category and consists, for 
example, of intrinsic error introduced by the biological 
components. In this regard, regulation at both the gene 

20 expression and protein activity levels can be 
substantially altered due to apparent negligible 
experimental differences in the treatment of a biological 
sample. A specific example is where gene expression 
changes due to the use of different batches of the same 

25 media during the course of an experiment. Such 

biological error produces measurable differences in the 
level of an analyte such as an expressed gene. 

Another category is the extrinsic error 
introduced through experimental manipulation. For 
30 example, differences in sample preparation, analyte or 
probe labeling efficiency, hybridization or binding 
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conditions, synthesis of probes, batches of solid-phase 
substrate and detection efficiency introduce variations 
in the determination of a measured analyte, even though 
all components and processes can be controlled so as to 
5 result in apparent negligible differences. Nevertheless, 
measurable differences in observed analyte signal occur 
due to the introduction of such error. 

Similarly, for non-biological settings the 
methods of the invention are applicable for determination 
10 of true signals from measured analytes in essentially any 
process or steps thereof for which a quantitative 
determination or comparison of a measurable component is 
desired. 

The above exemplary, and other forms of error 

15 all affect the perceived amount of a measured analyte 

through the introduction of fluctuations in the observed 
signal. Assessing the true signal of the analyte, 
independent of such fluctuations, allows direct 
comparison of analyte levels. Moreover, because the true 

20 signal of an analyte measurement can be determined, the 
methods of the invention provide a means for a direct or 
standardized comparison of analyte measurements both 
within an experimental system and between different 
systems. Given the teachings and guidance provided 

25 herein, essentially any analysis format known in the art 
can be used for such subsequent comparison of analytes 
once the true or mean signals are obtained. Therefore, 
the methods of the invention can be used to accurately 
and reproducibility determine the true signal of 

30 essentially any measurable analyte as well as used for 
the initial step in, for example, a comparative analysis 
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of the same analyte under different conditions, the same 
analyte under repetitive conditions or different analytes 
under the same conditions. 

As will be described further below, it is 
5 understood that the methods of the invention are equally 
applicable to both large and small sets of analyte 
samples and sets of measurements. Determination of the 
true signal for an individual sample is performed 
similarly as that for the determination of many, and even 

10 hundreds or thousands of samples. Similarly, the 

comparison of true signals for determination of relative 
amounts of an analyte between samples also is performed 
for two samples as it is for comparison of many sample 
pairs or higher order sets of multiple comparisons. 

15 Therefore, given the teachings and guidance provided 
herein, the number of true signals that can be 
simultaneous determined, or sets of samples that can be 
simultaneously compared for relative amounts of true 
signal is only limited by the available computational 

20 power. 

The methods of the invention for determining 
the true signal of an analyte can be applied to a variety 
of situations. For example, repeated measurements of the 
observed signal such as intensity x for one or more 

25 analytes can be obtained and subsequently used in the 
method of the invention to characterize the error and 
determine the significance value for each observed 
signal. For example, repeated observations of the signal 
associated with a single analyte such as, for example, 

30 the observed intensity of a single gene in a microarray, 
can be utilized in the methods of the invention to 
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monitor, for example, the variation introduced by two or 
more distinct conditions, the total error introduced over 
a given time or sporadic error introduced by any means 
including variation caused at any step in the protocol. 

5 The method of the invention provides a 

description of the relationship between an observed 
signal and a mean signal. The relationship specifies 
that the observed signal can be described as containing 
both an additive error term and a multiplicative error 

10 term. The error terms are a measure of variation in the 
observed signal. Parameters of the additive error term 
and the multiplicative error term set forth the 
characteristics or features of the error terms. These 
parameters are derived from statistical relationships 

15 well known in the art. Therefore, the error terms, and 
the parameters defining them, specify the noise of the 
analyzed system. Knowing the components and relationship 
of the noise with reference to the mean or true signal 
allows determination of the true signal given an 

20 empirically measured signal. 

The inclusion of both an additive error term 
and a multiplicative error term in the described 
relationship permits distinction of the true signal from 
the noise at a wide range of observed signals. For 

25 example, with a high observed signal, or high observed 
signal relative to the noise, the system noise can be 
primarily described by the multiplicative error term. 
Therefore, the true signal can be accurately 
distinguished from the noise by employing only a 

30 multiplicative error term in the method of the invention. 
In contrast, where the observed signal is low, or low 
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relative to the noise, the influence of the additive 
error in describing the noise becomes substantially more 
prominent. Maintaining this error term in the described 
relationship at low observed signals enhances the 
5 accuracy in distinguishing the true signal from the 
noise. Similarly, at intermediate observed signal 
ranges, both the additive and multiplicative error terms 
substantially influence the description of the noise and 
inclusion of both will yield enhanced results in 

10 distinguishing the true signal from the noise using the 
described relationship in the method of the invention. 
Therefore, including both the additive and multiplicative 
error terms in the description of the relationship 
between the observed signal and the true signal results 

15 in more accurate and predictable performance of the 

method of the invention at all ranges of observed signal. 

However, utilization of both the additive and 
multiplicative error terms in the methods of the 
invention is not always necessary. As described above, 

20 if the user knows or can determine that the observed 
signal is high relative to the limits of detection or 
relative to the noise, then determination of the true 
signal can be accurately made by inclusion of only the 
multiplicative error term. In such circumstances, the 

25 additive variation will be small or negligible compared 
to the observed signal and is included in the described 
relationship as an example where one or more of the error 
term parameters, such as the standard deviation of the 
additive error term, is set to zero. Similarly, where 

30 the signal is low but the variation is also known, or can 
be determined to be small, in like manner the additive 
error term also can be omitted without substantial affect 
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on determination of the true signal. Determination of 
the true signal also can be accurately made by inclusion 
of only the additive error term. For example, applying 
only the additive error term in the described 
5 relationship can be useful for measuring the error in the 
variation of the background of a system. Given the 
teachings and guidance described herein, those skilled in 
the art will know, or can determine, whether 
determination of a true signal can be made, or is 
10 desirable, utilizing both the additive and multiplicative 
terms in the described relationship employed in the 
method of the invention. 

For each analyte, the method of the invention 
provides a relationship between the observed signal and 
15 the mean signal which can be described as follows: 

where each measurement j equals 1 through M and each 
analyte i equals 1 through N, and where x^ is the 
observed signal and u xi is the mean or true signal. For 
each analyte and measurement, the multiplicative error 
term, £ x , and the additive error term, S x , can be 
obtained, for example, from a normal distribution with 
mean zero and standard deviation a ex and a 5x , 
respectively. One advantage of the above described 
relationship is that the multiplicative and additive 
errors can be independent of one another. Additionally, 
the additive and multiplicative error terms can be 
derived from a variety of univariate distributions, 
including, for example, a parametric distribution, a 
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univariate normal distribution, a t-distribution or a 
gamma distribution . 

For determining the true signal of an analyte, 
where the observed signal x^ is described by a univariate 
5 distribution with the parameters \i xl and o xl , the error 
model specifies two analyte-independent parameters, which 
together consist of the system parameter (3), and a mean 
signal |i x for the analyte. The system parameter |3 
describes the noise in the observed signal and consists 

10 of the above described standard deviation of the 

multiplicative error with respect to the mean (a sx ) and 
the standard deviation of the additive error with respect 
to the mean (a 5x ) . A particular feature of the above 
relationship, or error model, is that it specifies both 

15 the mean signal and noise such that the estimate of the 
signal describes the structural features of the noise. 
Therefore, the system parameter specifies the properties 
of both the multiplicative and additive error and can be 
independent of the mean signal. Moreover, the error 

20 terms specified in the model can be independent of one 
another . 

Modifications can be incorporated into the 
general description of the relationship between the 
observed signal and the mean signal set forth above and 

25 below which do not alter the relationship of the additive 
or multiplicative error terms with respect to the true 
signal or their properties in specifying the structure of 
the noise. Such modifications are exemplified with 
reference to the description specifying the relationship 

30 between the observed signal and the true signal set forth 
above, but are similarly applicable to the description 
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specifying the relationship between observed and true 
signals for comparison of two or more signals. The 
modifications can include, for example, inclusion of 
functions, augmentations or addition of terms, 
5 simplification or removal of terms and transformation of 
variables. Depending on the origin of the signal data or 
the desired use, one or more of such modifications can be 
employed to generate alternative forms of the described 
relationship appropriate for application to a wide 
10 variety of data sets. These modifications as well as 

others are well known to those skilled in the art and are 
applicable in the method of the invention. 

For example, the description specifying the 
relationship between the observed and true signal can be 

15 modified by inclusion of a function such as f(iu xi )£ xij 
where f is a function that describes how the mean 
sensitive component of variability varies as the mean 
varies. The function can do so simply by multiplying the 
mean signal by £ xij , or it can do so by multiplying e Ki3 by 

20 other terms related to the mean, in addition to the mean 
or together with the mean. Additionally, the system 
parameters also could be chosen as a function of the mean 
parameter }i xi . For example, and with respect to the 
expanded relationship set forth below describing the 

25 comparison of two or more true signals, the system 

parameters p £ and p 6 can be chosen as a function of the 
mean parameters p. xi and p. yi . With either of the above 
exemplarily functions, the system parameters would change 
according to principles well know to those skilled in the 

30 art to reflect the joint properties of the error of the 
system given the teachings and guidance provided herein. 
For example, the function "f" can be chosen to be a 
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polynomial of low order and the system parameter would be 
enlarged to include the coefficients of these enlarged 
polynomials . 

The description specifying the relationship 
5 between the observed and true signal also can be modified 
by augmentation. For example, terms can be added which 
include constants, second order or even higher order 
terms which do not alter the relationship of the additive 
or multiplicative error terms with respect to the true 

10 signal or their properties in specifying the structure of 
the noise. A specific example of the addition of a 
constant is Xij = u xi + u xi £ xij + 5 xij + C, where C is a 
global parameter which allows, for example, translation 
of the relationship along selected axes. Shifting the 

15 distribution by a constant can be useful, for example, in 
the normalization process to better fit the data as a 
whole. Additionally, a specific example of the addition 
of a second order term is x i;j = u xi + (u xi + au 2 xi ) £ xij + 5 xij . 
A specific example of the addition of a higher order term 

20 is x i3 = p xi +(U xi + au 2 xi + Pu 3 xi ) s xij + 5 xij -. These latter 
two descriptions allow for curvature in the relationship 
between the mean signal and the standard deviation at 
medium-to-large signal intensities. 

Simplification or removal of terms has been 
25 described above, such as when there is a negligible 

amount of error. Removal of the corresponding error term 
can increase the accuracy of determining the remaining 
parameters and therefore the accuracy of determining the 
true signal. A specific example of a simplification 
30 modification where the additive error has been removed is 

Xij = U x i + Uxi^xij* 
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Transformation of variables is yet another 
modification which can be performed that does not alter 
the relationship of the additive or multiplicative error 
terms with respect to the true signal or their properties 
5 in specifying the structure of the noise. For example, 
because some signal measurements can be distributed over 
a large range of values, including many orders of 
magnitude, it can be useful to transform the raw signal 
measurements into logarithms. For this transf ormation, 

10 the variables x ijr or for example y i3 in the relationship 
set forth below, can be redefined in terms of other 
variables such as s and t. Specifically, define s = 
log(x i j) and take the log of both sides of the equation: 
log(x ij ) = log(x i;j - u xi + U si £ xij + 5 xij ) . In the specific 

15 case where the additive error is small, the above 

equation reduces to: log{x i:j ) = log(p. xi ) + log (l+£ xij ) . 
Substituting s = log(x i;j ), this equation relates the 
sample value of s to the mean of s plus some additive 
error f = log(l+£ xij ), as in: s = ji s + f. Other 

20 transformations include, for example, exponentiation (s = 
£ n xij ) or polynomial transformations (s = ax n ij ) . 

The methods of the invention employ the above 
error model to determine, by formal estimation, the mean 
signal of an analyte from a set of measurements of an 

25 observed signal by using a maximum likelihood approach. 
To estimate the mean signal, the observed signal should 
be measured at least twice ( j=2) , obtaining two separate 
values and allowing for a more accurate computation of 
the system parameter and mean signal. However, a larger 

30 number of analyte measurements, where j is greater than 
2, results in further refinements of true signal 
determination. For example, as shown in Example I, 
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increasing the number of measurements from two to four 
per analyte results in beneficial enhancements in true 
signal determination. Therefore, the number of 
measurements of a particular analyte can be a few or many 
5 times, including for example, about 2, 3, 4, 5, 10, 20, 
50, 100 or more sample measurements. Although as few as 
two measurements is sufficient to accurately determine 
the true signal of an analyte, the actual number of 
measurements will vary depending on the need and 

10 confidence requirement of the user. For example, the 

confidence in true signal determination can be increased 
in analyte samples exhibiting inherently greater 
variation by compensating for the greater experimental 
error through increasing the number of sample 

15 measurements. Sample measurements can be derived, for 

example, from independent samples, replicates of the same 
sample that are independently measured, repeated 
measurements of the same sample or any combination 
thereof. 

20 Once the signal has been measured for one or 

more analytes, the observed signals can be subjected to a 
variety of statistical methods well known in the art to 
prepare the raw data for maximum likelihood analysis. 
Such methods include, for example, standardization and 

25 filtering techniques. Briefly, non-specific background 
can be subtracted to produce, for example, the observed 
signal x' . Moreover, depending on the need, the data 
measurements can be, for example, normalized to have 
comparable medians and extreme signals within a set of 

30 multiple measurements that are artif actually outside the 
signal range of its partners can be removed. Such 
modified values for the observed signal are similarly 
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applicable in the methods of the invention for 
determining the true signal of an analyte. Therefore , 
the error model of the invention additionally accounts 
for the influence of multiplicative and additive errors 
5 on the observed signals and provides a relationship 

between an observed signal x' , and the corresponding mean 
or true signal. 

As will be described further below in context 
of a comparing relative differences between two or more 

10 true signals, once obtained for any particular set of 
analyte measurements, the observed signal x or x' is 
analyzed by, for example, maximum likelihood probability 
for determination of its mean signal. In addition to a 
maximum likelihood approach, other approaches are known 

15 in the art to determine, by formal estimation, the mean 
signal from a set of observed measurements, including, 
for example, Quasi-Maximum Likelihood and Generalized 
Method of Moments. 

In addition to determining the true signal of 
20 an analyte, the methods of the invention also can be 
utilized to determine relative amounts of an analyte 
between samples. Briefly, following the methods 
described above for determination of a true signal for an 
individual analyte, for comparison of relative amounts of 
25 two or more analytes, observed signals are measured for 
each analyte and the corresponding true signals 
determined by probability likelihood analysis. The 
resultant true signals are then formally assessed by, for 
example, a difference indicator to determine relative 
30 levels. In this embodiment, for example, the methods of 
the invention identify true signals that are 
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significantly unequal, thus representing different 
amounts of analytes between the compared samples. 

The methods of the invention allow relative 
comparison of true signals between two analytes or pairs 
5 as well as between multiple analytes or sets. As 

described previously, the analytes to be compared can be 
can be different, or they can be substantially the same 
species of analyte but subjected to distinct conditions 
or obtained from distinct sources. Briefly, samples 

10 harboring analytes to be compared are referred to herein 
as sample pairs or sets. True signals resulting from 
each observed analyte signal for a particular comparison 
are similarly referred to as mean signal pairs or mean 
signal sets. Similarly, the true signals being compared 

15 for substantially the same analyte species derived from 
different conditions or sources is referred to herein as 
mean signal pairs per analyte and mean signal sets per 
analyte . 

By reference to comparison of two analytes, for 
20 the determination of relative amounts of an analyte 
between samples the observed signal and mean signal 
within a sample pair can be described by the following 
relationship: 

Xij = Uxi + Uxi^xij + 5 X ij/ and 
25 y Aj = u yi + u yi £: yij + 5 yij 

where each measurement j equals 1 through M and each 
analyte i equals 1 through N; where x ±j and y Xj are the 
observed signals, and where |i xi and \x yi are the mean 
signals. For each pair of analytes and measurements, the 
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multiplicative error terms, £ xij and £ yij , can be obtained, 
for example, from a bivariate normal distribution with 
mean zero and standard deviations a EX and a £y , and 
correlation p e . Similarly, the additive error terms, 5 xij 
5 and 5 yij also are drawn from a bivariate normal 

distribution with mean zero and standard deviations a 5x 
and a 5y , and correlation p 5 . Aside from the correlations 
already described, the error terms for a particular 
analyte i can be independent, that is, the multiplicative 

10 error terms (£ xi and e yl ) can be independent of the 

additive error terms (5 xi and S yi ) , and the error terms 
for analyte i (£ xi , e yl , 5 xi , 5 yi ) can be independent of the 
error terms for analyte j(£ xj , £ yj , 5 xj , 5 yj ) when j does 
not equal i (j*i). Additionally, the additive and 

15 multiplicative error terms can be derived from a variety 
of other bivariate distributions, including for example, 
a parametric distribution, a bivariate normal 
distribution, a t-distribution or a gamma-distribution, 
and, further, the independence assumptions can be dropped 

20 by including additional correlations in the system 
parameter (3 . 

The above described relationship between 
observed and mean signals for two analytes substantially 
parallels that described previously for an individual 

25 analyte. Therefore, this error model similarly provides 
the advantage of allowing multiplicative and additive 
errors to be independent of one another. Similarly, the 
above described error model can be applied by analogy to 
determination true signals for multiple analytes, 

30 including three or more analytes. For example, similar 
mean signal, multiplicative and additive error terms for 
analyte z can be described in a third equation. 



Additionally, higher order comparisons and error models 
can additionally be described using the teachings and 
guidance provided herein. 

For determining the true signal of an analyte 
pair, where, for example, the observed signals x Xj and y ±j 
are described by a bivariate distribution with the 
parameters u xi , u yi , a xi , a yi and p xiyi the error model 
specifies six analyte-independent parameters, which 
together consist of the system parameter (3=(a £x , a ey , p e , 
a 5x , a 5y , p 5 ) , and a mean signal pair, (u xi , u yi ) for the 
analyte. As with the univariate distribution described 
previously, the system parameter p for the bivariate 
distribution similarly describes the noise in the 
observed signal and consists of the above described 
standard deviation and correlations. Briefly, the 
analyte-independent parameters of the system include the 
standard deviation of the multiplicative error with 
respect to the mean of signal x (a ex ) , the standard 
deviation of the multiplicative error with respect to the 
0 mean of signal y (a ey ) , a correlation of the 

multiplicative error for the mean of signals x and y 
(p e ), the standard deviation of the additive error with 
respect to the mean of signal x (a Bx ) , the standard 
deviation of the additive error with respect to the mean 
5 of signal y (a 6y ) and a correlation of the additive error 
for the mean of signals x and y (p 5 ) . As described 
previously, one particular feature of the above 
relationship, and with the error models of the invention, 
is that it specifies both the mean signal and noise such 
■0 that the estimate of the signal describes the structural 
features of the noise. Therefore, the system parameter 
specifies the properties of both the multiplicative and 



additive error and can be independent of the mean signal. 
Moreover, the error terms specified in the model can be 
independent of one another. 

To determine, by formal estimation, the mean 
signal pairs of a sample pair, the observed signals x and 
y should be measured at least twice as described 
previously. Once the signals have been measured for 
analytes within one or more sample pairs, the raw data 
can be prepared for maximum likelihood analysis to 
produce, for example, two signals x' and y' . For 
analysis of more than two analytes within a sample pair, 
standardization and filtering methods can similarly be 
used to produce, for example, signals z' and the like for 
sample sets. These methods and others well known in the 
art for processing raw data into useful statistical form 
are particularly appropriate when analyzing multiple 
observed signals of sample pairs and sets in order to 
provide meaningful comparisons by, for example, 
normalization of divergent scales for the initially 
measured signals. Such modified values for the observed 
signals are similarly applicable in the methods of the 
invention for determining mean signal pairs, mean signal 
pairs of an analyte and mean signal sets. Therefore, the 
error model of the invention additionally accounts for 
the influence of multiplicative and additive errors on 
the observed signals and provides a relationship between 
observed signals x' , y' , z' and higher numbers of like 
comparisons, and the corresponding true signals. 

For any of the error models described above, 
once an observed signal, observed signals within a sample 
pair or sample set are obtained, the mean signal (u) and 
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the system parameter (3) can be determined or selected 
by, for example, a non-linear optimization algorithm. 
Such statistical optimization procedures are well known 
in the art and can be applied to, for example, individual 
5 observed analyte signals, observed signals for a single 
sample pair and to observed signals for two or more, 
including, for example, hundreds, thousands or ten 
thousand or more signals for sample pairs or sets. The 
number of optimizations that can be performed is 
10 coextensive with the number of analyte signals or higher 
order sets that can be measured and the computing power 
available in the art. 

Similarly, and in addition to non-linear 
optimization algorithms, any general optimization 

15 procedure for non-linear equations can be used to 

determine or select the mean signal pair (u) and a system 
parameter ((3) for each sample pair including, for 
example, Gradient Descent, Newton-Raphson and Simulated 
Annealing. For example, The Gradient Descent method is 

20 based upon selecting, at each iterative step, the 
direction in multidimensional space for which the 
objective function initially changes at the fastest rate, 
and subsequently choosing an appropriate distance to move 
in this direction at that iterative step. The Newton- 

25 Raphson method is based on a linear approximation to the 
first-order conditions, which may be numerically 
estimated, that set to zero the partial derivatives of 
the objective function with respect to the parameters 
being estimated. The Simulated Annealing method is based 

30 upon making random changes, which become smaller 
throughout the iterations, in the parameters being 
estimated and subsequently deciding probalistically 
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whether or not to keep these changes, thereby seeking an 
optimum while maintaining the ability to escape from a 
suboptimal local optimum in^ order to seek a better 
solution. 

5 Further, the methods of the invention also 

allow the mean signal and system parameter to be provided 
based on previously determined or estimated values rather 
than calculated de novo. For example, in routine or 
familiar procedures, the user can have prior knowledge of 

10 beneficial or optimal estimates that can be used to 

calculate enhanced values for the probability likelihood 
or which more efficient convergence to a maximum 
probability likelihood. Therefore, the mean signal pair, 
including the mean signal pair per analyte, for example, 

15 (ia) and a system parameter ((3) for each sample pair can 
be determined or provided and then subsequently compared. 
As will be described further below, comparison of mean 
signals, mean signal pairs and higher order sets can be 
performed, for example, by identification of 

20 significantly unequal mean signals using well known 
methods in the art such as statistical difference 
indicators . 

In one embodiment, the mean signal and system 
parameters are estimated using maximum likelihood 

25 estimation. The maximum likelihood function provides, 
for example, a framework for the formal estimation 
process, while recognizing the structure of the random 
noise in the system. By modeling patterns of randomness, 
the maximum likelihood estimation process can better 

30 separate and estimate the signal. The method of the 

invention provides likelihood functions using estimates 
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for the true parameters by utilizing standard 
optimization procedures as described herein. One 
advantage of the methods of the invention is that, if 
desired, the error terms can be independent of one 
5 another. Moreover, each mean signal within a mean signal 
pair or set also can be independent with respect to each 
other. These characteristics allow for the independent 
optimization of the system parameter and mean signal. 
Therefore, the efficiency of optimization can be 
10 significantly increased for a large number of analytes, 
for example, through the optimization of the system 
parameter and mean signals in subsets. 

Briefly, observed values are measured and, 
subsequently, the system parameter (3) can be selected to 

15 enhance the probability likelihood given the observed 
signal. Similarly, for each analyte, mean signal pairs 
can be selected to enhance the probability likelihood 
given the system parameter (p) . The mean signal pair and 
system parameter can be determined at the same time, or 

20 alternatively, the mean signal can be determined prior to 
the system parameter and then subsequently used to 
determine the system parameter. Conversely, the system 
parameter can be determined prior to the mean signal and 
then subsequently used to determine the mean signal. As 

25 described further in Example I, this procedure can be 
reiterated one or more times until the mean signal pair 
per analyte (u) and a system parameter (3) converge. 
With each selection of values and reiteration of the 
optimization procedure, the calculated mean is enhanced 

30 in the direction of the true signal for that analyte, 
pair or set. In addition to maximum likelihood 
estimation, probability likelihood values for system 
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parameters and mean signal can be estimated using other 
modeling techniques known in the art including, for 
example, Quasi-Maximum Likelihood and Generalized Method 
of Moments. 

5 For comparison of the relative levels of two or 

more true signals, after the system parameter and mean 
signal have been determined, the methods of the invention 
provide for identification of mean signal pairs that are 
significantly unequal, representing different amounts of 

10 analytes between the compared samples. The error models 
and methods of the invention take into account the 
observation that x and y variances and x-y correlation 
increase with increasing values of x and y. Based on 
these empirical observations, the methods of the 

15 invention utilize a likelihood ratio test to identify 
analytes whose true signals ]i x and y y are unequal. For 
example, in the case of RNA expression analysis, analytes 
with unequal mean signals have different copy numbers of 
the measured mRNA analyte in the two cell populations 

20 under comparison, or in other words, are 

differentially-expressed. Such methods for assessing 
significantly unequal mean signals are well known in the 
art and are described further below in the Examples. 
Thus, the methods of the invention provide a difference 

25 indicator for comparison of true signals and therefore 
relative amounts of two or more analytes. Additionally, 
when used in combination with known analyte standards, 
the methods of the invention can be employed to 
quantitate the amount of a test analyte by comparison of 

30 its true signal with that of one or more known standards. 
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The methods of the invention can be utilized 
for determining the true signal of an analyte or for 
comparing the relative levels of two or more true analyte 
signals in a variety of different formats and modified 
5 procedures. For example, observed signals for one or 
more analytes, sample pairs or sample sets can be 
measured independently, such as in series, or 
simultaneously, such as in parallel. Moreover, different 
observed signals can be measured, for example, from 

10 independent samples, the same sample or from independent 
samples that have been pooled to reduce the total number 
of samples which are to be manipulated. The number of 
different observed signals which can be measured from a 
single sample or pooled sample will depend, for example, 

15 on the number of unique detection labels which can be 
employed to uniquely measure each different analyte 
within the sample. Corresponding mean signals, mean 
pairs or mean sets can similarly be determined from the 
observed signals in series or parallel, for example. 

20 Additionally, the measurements of observed signals and 
determination of mean signals can be multiplexed with 
ongoing measurements and determinations proceeding 
simultaneously in series or parallel, such as in an 
automated system, for example. 

25 Various modification can be made to the 

procedure described above for determining or comparing 
true signals which enhance the description of the noise 
and therefore, further increase the accuracy of 
distinguishing the true signal from the noise. For 

30 example, variation of a reference signal can be captured 
or incorporated into the analysis. In this specific 
example, two or more observed signals to be compared are 



first independently compared to a reference signal to 
determine, for example, the system parameters or mean 
signal pairs for each test-reference comparison, A 
probability likelihood can then be generated from the 
product of the terms for each initial test-reference 
comparison, to describe, for example, p L and (3 2 . These 
system parameters obtained with respect to the test- 
reference comparison can then be used, for example, to 
determine the mean signal pairs or sets for the two or 
more observed signals to be compared. Briefly, and as 
described further below, a likelihood is then established 
as the product of Mf^, u 1 ^, ]i\ ± ) and Mf^, u 2 xi , U 2 yi ) . A 
statistical difference indicator can then be applied, for 
example, constraining u 1 ^ and p 2 xi , as well as ]i 1 yi and u 2 yi 
to be equal or not equal to each other as described 
previously. For the specific example where y represents 
the reference sample, then u x yi and |i 2 yi would be 
constrained to be equal. Variation can be captured from 
one or more reference signals alone or in combination. 
0 Additionally, using the teachings and guidance provided 
herein, other methods well known to those skilled in the 
art which enhance the description of the signal or noise 
can additionally be incorporated into, or used in 
conjunction with the methods of the invention. 

5 Therefore, the invention provides a method of 

determining relative amounts of an analyte between 
samples. The method consists of: (a) obtaining a 
reference signal; (b) obtaining observed signals x and y 
for an analyte within two or more sample pairs; (c) 

0 determining system parameters (3 2 ) for a sample pair 

comprising said observed signals x or y and said 
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reference signal that provide a probability likelihood of 
said occurrence given said observed and reference 
signals, said observed and reference signals being 
related to said mean signal by an additive error (5) and 
5 a multiplicative error (£), where said system parameter 
specifies the properties of said additive error and said 
multiplicative error; (d) determining mean signal pairs 
(]i lf \x 2 ) for said sample pair comprising maximizing a 
product of terms for said probability likelihood of said 
10 sample pair of observed signals x or y and said reference 
signal for said analyte, and (e) selecting a mean signal 
]i x or u y that provides a maximum probability likelihood of 
occurrence given said observed signals and system 
parameters 3 X and |3 2 . 

15 The invention also provides a method of 

determining relative amounts of large numbers of analytes 
between samples. The method consists of: (a) obtaining 
observed signals x and y for a plurality of immobilized 
analytes within two or more sample pairs; (b) determining 

20 a mean signal pair per analyte and a system parameter 

(3) for each sample pair that provides a maximum 
probability likelihood of occurrence given the observed 
signals, the observed signals being related to the mean 
signal by an additive error (5) and a multiplicative 

25 error (£), where the system parameter specifies the 

properties of the additive error and the multiplicative 
error, and (c) identifying one or more mean signal pairs 
per analyte that is significantly unequal. The method is 
applicable, for example, to nucleic acid and polypeptide 

30 analytes using immobolized array formats. 
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The methods of the invention are applicable for 
determination or comparison of true signals in a wide 
variety of systems. Various detection methods for 
numerous analytes are well known to those skilled in the 
5 art. All that is needed to practice the methods of the 
invention are measurable quantities of an analyte in a 
data form that can be calculated as a mean. 

In biological systems , for example, detection 
of a nucleic acid analyte can be by any of a variety of 

10 detection methods well known to those skilled in the art. 
Such methods include, for example, gels, blots, 
capillaries and microarray formats. In addition to 
nucleic acid microarrays or chips, the methods of the 
invention further can be applied to determine the true 

15 signal of polypeptide spotted on a chip. The 

construction of glass chips or other substrates spotted 
either with chemicals to bind polypeptides or with known 
antibodies can be constructed and the bound polypeptide 
analyte can be detected, for example, by a mass 

20 spectrometer. Moreover, detection of a polypeptide 
analyte also can be by any other of a variety of 
detection methods well known in the art, including, for 
example, gels, blots, capillary and FACS formats. In 
addition, analytes other than nucleic acids and 

25 polypeptides can be detected by methods known in the art 
such as spectroscopy and laser-assisted techniques. The 
detection method and, consequently, the visualization 
technique that yields the observed signal will depend on 
a variety of factors such as the nature, amount, 

30 stability and purity of the analyte. 
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Microarray hybridization and fluorescent 
detection is one well known method for analysis of large 
numbers of nucleic acid analytes. Currently, arrays with 
more than 250,000 different oligonucleotide probes or 
5 10,000 different cDNAs per square centimeter can be 

produce in significant numbers. Although it is possible 
to synthesize or deposit DNA fragments of unknown 
sequence, generally, microarray-based formats utilize 
specific sequences attached to a solid substrate such as 
10 glass, plastic, silicon, gold, a gel or membrane, beads, 
or beads at the ends of fibre-optic bundles. Such 
formats allow for parallel hybridization and simultaneous 
detection of a large number of indexed, surface-bound 
nucleic acid probes. 

Nucleic acid arrays are generally produced by 
either robotic deposition of nucleic acids such as PCR 
products, plasmids or oligonucleotides, onto a glass 
slide or in situ synthesis using, for example, 
photolithography of oligonucleotides. After 
hybridization of labelled samples to the spotted or 
synthesized probes, the arrays are scanned and a 
quantitative fluorescence image along with the known 
identity of the probes is used to detect the presence of 
a particular molecule above thresholds based on 
background and noise levels. 

Various methods for preparing labelled material 
for measurements of gene expression microarrays are well 
known in the art. For example, the RNA can be labelled 
directly, using a psoralen-biotin derivative or by 
30 ligation to an RNA molecule carrying biotin, labelled 
nucleotides can be incorporated in cDNA during or after 
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reverse transcription of polyadenylated RNA; or cDNA can 
be generated that carries a T7 promoter at its 5 T end. 
In the last case, the double-stranded cDNA serves as 
template for a reverse transcription reaction in which 
5 labelled nucleotides are incorporated into cRNA. 
Commonly used labels include the fluorophores 
fluorescein, Cy3 or Cy5, or nonf luorescent biotin, which 
is subsequently labelled by staining with a fluorescent 
streptavidin conjugate. Generally, cDNA from two 
10 different conditions is labelled with two different 

fluorescent dyes such as Cy3 and Cy5, and the two samples 
are co-hybridized to an array. After washing, the array 
is scanned at two different wavelengths to detect the 
relative transcript abundance for each condition. 

Another quantitation method which is useful for 
determining expression levels of polypeptide analytes is 
the isotope-coded affinity tag (ICAT) method (Gygi et 
al., Nature Biotechnol. 17:994-999 (1999)). 
Specifically, ICAT involves labeling two analyte samples 
differently by using stable isotopes, loading them into a 
mass spectrometer, and measuring the ratio of the two 
labels and thus the relative mass. ICAT can make any 
separation method, including HPLC and capillary 
electrophoresis, quantitative and, rather than using a 
ratio-based comparison, the methods of the invention can 
be applied to any of these separation methods to 
determine the true signal of a polypeptide analyte or the 
relative amounts of an analyte between samples. 

Additionally, measurement of an analyte signal 
30 can be by a variety of other methods well known in the 
art, including, for example, light emission, 
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radioisotopes, and color development. Briefly, detection 
can involve methods such as radioactive labeling of the 
analyte using metabolic labeling in an appropriate cell 
or in vitro labeling by RNA transcription or by coupled 
in vitro transcription-translation with appropriate 
radioactive amino acids. Additionally, covalent 
modification with a radioactive or fluorescent substrate 
using an appropriate enzyme or chemical modification can 
be employed. Moreover, an analyte can be covalent ly 
modified by incorporating a chemical moiety capable of 
being detected. For example, green fluorescent protein, 
Cy3, Cy5 and other fluorophores can be covalently 
attached to a polypeptide analyte. Similarly, biotin can 
be covalently attached to a polypeptide analyte and 
subsequently detected by streptavidin using detection 
methods known in the art. Other methods also can involve 
fusion of an appropriate detection molecule to the 
analyte. For example, the analyte can be fused to 
luciferase and detected by light emission or can be fused 
to lacZ and detected by appropriate calorimetric 
detection. 

The methods of the invention have utility for a 
variety of applications. Although a standard microarray 
compares only two populations, a greater number can be 
cross-compared by hybridizing labeled probe, such as cDNA 
prepared from each cell population of interest, to that 
of a common reference population. The methods of the 
invention can thus be used to determine genes 
differentially-expressed between any two populations, 
even if they have not been directly involved together in 
a single hybridization experiment. 
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The error model of the invention does not 
distinguish between repeated samples drawn from multiple 
spots on a single array versus repeated samples drawn 
from multiple hybridizations to different arrays. 
5 Because multiple spots within an array show less 
variability and more dye-to-dye correlation than do 
multiple spots observed over several arrays, the error 
model of the invention can be applied to distinguish 
between these two types of sampling, resulting in a more 

10 sensitive or accurate likelihood ratio test. Systems 
which involve more than one level of sampling are well 
known in the art and can be addressed by utilizing a 
nested design model as described by Dunn and Clark, 
Applied Statistics: Analysis of Variance and Regression 

15 (John Wiley & Sons, Inc., New York, 1987), which is 
incorporated herein by reference. 

The methods of the invention further can be 
utilized to place a confidence interval on the true 
signal difference between two analytes. In this 
20 embodiment, rather than testing the hypothesis that U x =U y , 
the range K(u x -u y )<h or the range 1< (u x /u y ) <h is 
determined for each analyte. 

In another embodiment, the methods of the 
invention can be utilized to quantify, compare, and 

25 ultimately reduce the error introduced by each stage of 
an array process. Therefore, the methods of the 
invention can be used for quality control in a large 
variety of processes and settings. For example, as shown 
in Example II, system parameters and mean signals can be 

30 compared for replicate spots on one array versus a single 
spot observed over multiple array hybridizations (see 
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also Table 2) . It is understood that this embodiment of 
the method of the invention can be expanded to quantify 
several different levels of variation, such as variation 
due to cell culture, RNA preparation, labeling, or 
5 hybridization. Moreover, it can be expanded to other 
biological assay systems as well as non-biological 
systems. Thus, the method of the invention can be 
utilized to identify sources of variation that contribute 
to the overall error of the system. 

10 The methods of the invention can be extended to 

a wide range of biological data involving comparisons 
between multiple measurements and can be advantageously 
utilized to determine differential gene expression based 
on studies with fluorescent or radioactive-labeled cDNA 

15 hybridized to gene clones spotted on membranes. 

Furthermore, the methods of the invention are applicable 
to large scale genotyping of human polymorphisms, where 
normal DNA is cut into small fragments, labeled, 
transferred onto a microchip and subsequently hybridized 

20 with labeled samples of normal and polymorphic DNA. 

Because the observed quantities of polypeptide expression 
per gene are analogous to fluorescent signals observed in 
a microarray experiment and are correlated, the methods 
of the invention can be practiced with technologies for 

25 comparing levels of polypeptide expression between two 
cell populations, for example (Gygi et al., Mol. Cell 
Biol. , 19:1720-1730 (1999), supra. Thus, the method of 
the invention can be advantageously utilized for 
describing measurements obtained in various technologies 

30 including those pertaining to, for example, genomics and 
proteomics . 
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For example, the method of the invention can be 
applied to proteomics where increased sensitivity of 
sequencing methods and mass spectrometry allow for 
determination polypeptide expression profiles. The 
5 methods of the invention can be advantageously used to 
determine relative amounts of polypeptide based on, for 
example, virtual 2-D profiles obtained by linking of 
isoelectric focusing gels with mass spectrometry. 

It is understood that the observed signal 
10 depends on the method of detection. For example, in the 
case of a microarray, the amount of hybridization can be 
quantified by, for example, optical imaging or laser 
scanning to observe the emitted light intensity. The 
observed signal also can be obtained by other 
15 visualization techniques based on the nature of the 

analyte as well as the assay and include, for example, 
chemi luminescence and fluorescence imaging systems, and 
mass spectrometry. These and other methods are well 
known in the art and can be employed for the detection of 
20 an observed signal in the methods of the invention. 

EXAMPLE I 

Development of an Error Model of the V ariability Observed 
over Repeated Observations of Intensit ies for Genes 
Represented on a DNA Microarray 

25 

This example describes development of a 
maximum-likelihood test for the variability observed over 
repeated observations of intensities for genes 
represented on a DNA microarray. 
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Preprocessing of Microarray Data 

The amount of hybridization to each spot is 
quantified by scanning the array with a laser and 
observing the intensity of light emitted. Observations 
5 are made separately for the two dyes, such that two 
intensities x and y are observed for each spot on the 
microarray. This process does not behave 
deterministically in practice, such that multiple spots 
corresponding to each gene i hybridized under identical 
10 conditions will result in a distribution of intensities 
and y^ (1 < i < N; 1 < j < M) , where N is the number 
of genes represented on the microarray and M is the 
number of spots observed for each gene. 

Spot intensities were extracted from a scanned 

15 image, then background-subtracted and normalized as 

follows: microarray images are processed with Dapple, a 
software tool developed for array spot finding and 
quantitation described by Buhler et al., Bioinf ormatics 
2000, which can be found at the URL: 

20 cs.washington.edu/homes/jbuhler/research/array, which is 
incorporated herein by reference. The Dapple software 
locates each spot and reports a separate median 
foreground intensity for each dye inside the spot area. 
The Dapple software also provides a local background 

25 intensity estimate for each spot and dye. The Dapple 

intensity estimates were subsequently smoothed by spatial 
filtering using a 7 spot by 7 spot median filter as 
described by Lim J.S. Two-Dimensional Signal and Image 
Processing (Englewood Cliffs, Prentice Hall, 1990) , which 

30 is incorporated herein by reference. Subsequently, the 
smoothed background was subtracted from the foreground of 
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each spot so as to produce the background-subtracted 
intensities x' and y' . 

In practice, X ! and y' have different scales 
and thus are not directly comparable- This situation can 
5 occur if the total amount of labeled cDNA is greater for 
one dye than it is for the other, if one dye incorporates 
more efficiently, or if the scanner has different 
sensitivities to the two dyes. Therefore, the 
intensities are normalized to have identical medians A 
10 within each array hybridization: 

x y ^ v J 

where X' denotes the median intensity of x' over all spots 
on a single microarray. If multiple array hybridizations 
are performed, normalization occurs independently for 

15 each and the resulting combined data set consists of data 
pairs {x 13 , y ±3 ) for gene i in repeat j. If three or more 
samples are available for a gene, these are filtered 
independently in x and y to remove outliers by Dixon's 
test with a=0.1 as described in Dunn and Clark, Applied 

20 Statistics: Analysis of Variance and Regression (2nd ed., 
Wiley and Sons, New York, New York, 1987), which is 
incorporated herein by reference. In addition, extremely 
high intensities outside the dynamic range of the array 
scanner in either color are removed. 

25 Formulation of the Error Model 

An error model summarizing the influence of 
multiplicative and additive errors on x and y has been 
formulated. In this regard, it has been consistently 
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observed that larger intensity measurements have a 
proportionately larger error over repeated samples. 

The data shown in Figure 1, which shows the 
increase of (A) standard deviation and (B) correlation 
5 with absolute level of intensity x' or y 1 , were obtained 
over 5 separate hybridizations with identically-prepared 
Cy3- and Cy5-labeled cDNA mixtures to test arrays 
containing 16 replicate spots per gene over 96 genes, 
resulting in a total of 80 samples for each of 96 genes. 

10 Figure 1 (C) shows the normal probability plot for the 80 
samples of x' pertaining to a single, representative 
gene. This plot is linear, indicating that these data 
are consistent with a normal distribution. The dotted 
line connects the 25th and 75th percentiles of the data 

15 and represents an approximate linear fit. 

As shown in Figure 1 (A) , larger intensity 
measurements have a constant coefficient of variation 
g^x' , as can be caused by variation in spot size or 
labeling efficiency from gene to gene. However, the 

20 variability does not tend to zero as x^O, likely due to 
variation in the measured background intensity. 
Furthermore, within genes, x and y are correlated and, in 
addition, larger intensities have a larger correlation, 
possibly due to errors introduced by spot-to-spot 

25 nonunif ormity or during the hybridization process which 
affect intensity measurements for both dyes 
simultaneously (see Figure IB) . Finally, as shown in 
Figure IB, samples of x and y for a given gene are at 
least approximately normally distributed, as assessed by 

30 a normal probability plot described by Dunn and Clark, 
supra r 1987. 
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Based on the observations described above, the 



background- subtracted, median-normalized intensities 
observed for each gene are related to their true (or 
mean) intensities by the following model: 



5 



= y*i + Uxi£*ij + 5 xijf and 

Yij = Pyi + Pyi^yij + 5 yij 



where (Ti xi ,U yi ) is the pair of true mean intensities for 
gene i. For each i and j , the multiplicative errors £ xij 

10 and £yij/ are drawn from a bivariate normal distribution 
with means 0, standard deviations g £x and G £y , and 
correlation p e . The additive errors 5 xij and 5 yij , are 
distributed analogously, with parameters a 5K , a 5y , and p 5 . 
Thus, multiplicative and additive errors are independent 

15 of one another but can each be highly correlated between 
x and y; in practice p E is large and p 5 is small. While 

and can be negative if the foreground is less than 
the estimated background for a spot, the true intensities 
TJL xi and ]i Y± must be non-negative. Consequently, the 

20 samples (x i5 and y^) are described by a bivariate normal 
probability density function p with parameters and ]i y± 
°x±f ° Y i and Pxi,yi/ where: 





p xi yj ~ 

(JxtCyi 



46 

The model depends on six gene-independent parameters 
3=(o £K , o ey , p c , G 6x , a 5y , p 5 ) and a mean pair per gene, 

U=[ (UxirUyi) / (Ux2,Vy2), • • • f (Uxn/ VW 1 for a total of 2N+6 
parameters. The probability density function for gene i 
5 is p=p(x ij7 y i:i |3, u xi , U yi ) . 

Parameter Estimation by Maximum Likelihood 

Since (3 and ]i are generally unknown, they can 
be estimated by using a maximum likelihood estimation 
(MLE) as described by Kendall and Stuart, The Advanced 
Theory of Statistics , Volume 2 (4 th ed., Macmillan 
Publishing Co., New York, N.Y., 1979), which is 
incorporated herein by reference. Likelihood functions, 
for gene i and over all genes, are respectively defined 
as : 

M 

M 

N 

^)=II Li{^flx h flyi) 
i=\ 

The MLE parameter values maximizing L, designated |3 and 
p., are estimates for the true parameters of the 
underlying statistical model. In general, these values 
can be found using standard optimization procedures as 
20 described by Press et al., Numerical Recipes in C: The 

Art of Scientific Computing (2 nd ed., Cambridge University 
Press, Cambridge, MA) . Because N can be large (3 and ]i r 
can be determined by optimizing subsets of parameters in 
separate stages: 
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(1) choose initial values for ]i r 

(2) select (3 to maximize L given current values of ]X f 

(3) for i=l, . N: select (u xi ,U yl ) to maximize L if given 
current values of (3, and 

5 

(4) repeat (2) and (3) until (3, y have converged. 

All stages of the optimization were performed using the 
procedure fmincon provided by Matlab and described by 
Coleman et al., Matlab Optimization Toolbox User's Guide 

10 (3 rd ed. f Mathworks, Inc., Natick, MA, 1999), which was 
incorporated herein by reference. The optimization was 
also implemented in C code, which produces comparable 
optimal parameters in substantially less execution time 
(less than 10 minutes on a Pentium III 500 for N=6000, 

15 M=4, as compared with 4-5 hours for the Matlab 

implementation) . In both cases, all parameters converged 
within 250 iterations of stages (2) and (3) and are 
insensitive to initial choices for 3 and ]i. 
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Significance Testing using Likelihood Ratios 

After the parameters have been determined for a 
given set of observations, it is of immediate interest to 
use the model to identify mean intensity pairs which are 
5 significantly unequal such that U X i*U yi , representing genes 
that are differentially expressed between the two cell 
populations. For each gene i, the generalized likelihood 
ratio test (GLRT) (Kendall and Stuart 1979) statistic X t 
is computed according to: 



10 



^ maxLi(/?,// ? //) ^ 



A, = -2\n 



M 

max Li(J},{ix,fi y ) 

V J 



Two maximizations are performed: in the 
numerator, the constraint ]i x =y. y =]a is imposed, while in the 
denominator the optimization is unconstrained. Under the 
null hypothesis that "|a x =ia y , (3 remains a consistent 
15 estimator when the constraint is imposed. 



In the case that U xi =Uyi, X L follows 
(asymptotically in M and N) a x 2 distribution with 
1 degree of freedom (DOF) , whereas if U xi *U yi , the value of 
\± is expected to be larger than would be obtained from 

20 random sampling of this distribution. To select 

differentially-expressed genes with a selection error of 
a, the false positive or Type-1 error rate, one would 
first determine the critical value A c , for which the x 2 
cumulative probability distribution is equal to 1-a, then 

25 select the set of all genes i for which X ± is in the 
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critical region A i >X c . The particular choice of a depends 
on the number of genes on the array and the selection 
error which the individual investigator is willing to 
tolerate. 



5 EXAMPLE II 

Identification of Genes Differentially-Expressed in 
Response to Galactose Stimulation of Yeast Cells 

This example describes application of the 
mathematical model of the variability observed over 
10 repeated observations of intensities for genes 

represented on a DNA microarray to the identification of 
genes differentially-expressed in response to galactose 
stimulation. 



Assembly of the Microarray 



15 In order to explore the performance of the test 

for differentially-expressed genes as shown in Example I, 
Saccharomyces cerevisiae cultures growing in the absence 
of galactose (YPR media) were compared to those growing 
in galactose-stimulating conditions (YPRG) using a DNA 

20 microarray of approximately 6200 nuclear yeast genes. 

The microarray was fabricated so as to consist of a large 
number of DNA spots on glass, each containing the full 
open-reading-frame sequence of a gene as reviewed by 
Lander, Nature Genetics 21: 3-4 (1999), which is 

25 incorporated herein by reference. 



Initially, mRNA contained in each of two 
populations of cells was extracted, reverse-transcribed 
into cDNA, and labeled with either Cy3 or Cy5 dye as 
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described below. Subsequently, the Cy3 and Cy5 dye 
preparations were combined and deposited on the 
microarray, where labeled molecules hybridize to the spot 
containing their complementary sequence. 

5 In order to obtain the mRNA to be 

reverse-transcribed into cDNA, wild-type yeast (BY4741) 
or a congenic gal80A strain were inoculated in 100 ml of 
either galactose-inducing YPRG media (1% yeast extract, 
2% peptone, 2% raffinose, 2% galactose) or non-inducing 

10 YPR media (1% yeast extract, 2% peptone, 2% raffinose) . 
Subsequently, cultures were grown at 30°C to a density of 
1-2 OD 600 , and total RNA was harvested by hot acidic 
phenol extraction as described by Ausubel et al., supra, 
(1995) . Poly-A purification from total RNA was performed 

15 using Ambion Poly (A) Pure mRNA Isolation Kits (Ambion, 
Austin, TX, catalogue #1915) . 

To assemble the DNA microarray a set of 
approximately 6200 known and predicted gene open reading 
frames from the yeast Saccharomyces cerevisiae (Research 

20 Genetics, Huntsville, AL) was amplified in separate IOOvlL 
PCR reactions in a 384-well plate format. The PCR 
conditions were optimized depending on the length of the 
template, but in general were as follows: Initially 95°C 
for 2 minutes; followed by 35 cycles of 94°C for 30 

25 seconds, 64°C for 30 seconds and 72°C for 2.5 minutes; 
and, finally, followed by 72°C for 5 minutes. The 
reaction products were subsequently purified over a 
Sephacryl S-500 spin column (Pharmacia, Uppsala, Sweden) . 
The purified product was then added to DMSO in a 1:1 

30 ratio. A Molecular Dynamics Generation III microarray 
robotic spotter was used to print the PCR products onto 



25mm by 75mm glass slides (Amersham, Piscataway, NJ, 
catalogue # RPK0328), which were subsequently spotted at 
50% humidity and immediately UV cross-linked at 50 mJ of 
energy. 

5 Complementary DNA synthesis and hybridization 

was accomplished as follows: 2iig anchored dT25 primers 
and 2]ag random 9-mer primers were added to 4p.g poly-A 
selected mRNA and allowed to anneal at 70°C for 5 minutes 
in a 12iiL volume. After 1 to 2 minutes on ice, 4pL 5x 

10 Superscript II buffer (Gibco) , 2pL 0.1M dTT, l]iL dNTP 
mix (lOmM dATP, dTTP, dGTP, and ImM dCTP) , ImM of either 
Cy3 or Cy5 fluorescent dye (Amersham, Piscataway, NJ) , 
and l]iL Superscript II reverse transcriptase were added. 
Reverse transcription occured at 42°C for 2 to 2.5 hrs in 

15 the dark. Subsequently, the RNA was hydrolyzed by 

heating at 94°C for 3 minutes, followed by addition of 
lp.L of 5M NaOH, and incubation at 37°C for 10 minutes. 
The pH was adjusted by the addition of IjjlL 5M HC1 and 5]iL 
1M Tris (pH 6.8) followed by cDNA purification through 

20 Millipore NAB plates (Millipore, Bedford, MA) . Dye 

incorporation was assessed by measuring absorbance at 550 
and 650 ran, and a sample aliquot containing about 40 pmol 
of dye is concentrated to less than 5 ]iL. Subsequent to 
labeling, purification, and concentration, Cy3 and Cy5 

25 samples were combined and suspended in 40 to 45 pL of 
hybridization solution containing 50% formamide, 5x 
Denhardt f s solution, 5x SSC and 0.1% SDS. The 
hybridization mixture was subsequently applied to the 
array slide beneath a coverslip and allowed to incubate 

30 in a sealed, humid chamber overnight for 16 to 18 hours 
at 42°C. The slide was then washed in 2x SSC/0.1% SDS 
for 5 minutes at 42°C, followed by a 5 minute wash in 
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O.lx SSC/O.1% SDS for 5 minutes at room temperature and, 
finally, two additional washes in 0.1 x SSC, each for two 
minutes. The slide was rinsed briefly in distilled water 
and immediately dried with compressed air* After 
5 hybridization and washing, the array slides were scanned 
using a scanning laser fluorescence microscope 
(Molecular Dynamics Generation III Scanner, Molecular 
Dynamics, Sunnyvale, CA) . 

Each gene was represented by two spots located 
on opposite sides of the array. A total of four (x,y) 
intensity pairs was obtained for each gene by performing 
replicate hybridizations to two of the above microarrays 
(N=6200, M=4), with x and y representing intensities in 
YPR and YPRG respectively. In the first hybridization, 
RNA from the YPR condition was labeled with Cy3 dye, 
while RNA from the YPRG condition was labeled with Cy5 
dye; in the second hybridization the reverse labeling 
scheme was used. The p and ]i values were determined for 
these data using our maximum likelihood approach, and the 
A x statistic was computed for each gene. Values for 3 
were as follows: 0.367, 0.391, 0.862, 89.6, 339.0, 0.319. 

In order to determine a reasonable choice for 
the critical value A c used to select 
25 differentially-expressed genes, a series of control 

experiments was performed in which two cell populations 
were cultured separately using identical strains and YPRG 
growth conditions. These two populations were compared 
as described before by obtaining a total of M=4 repeat 
30 samples per gene and determining values of (3, ]i and A. 
In general, these control data had fewer large values of 
A than did the YPR versus YPRG data, and followed a X 2 
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distribution as determined by a q-q plot. However, both 
data sets had significantly larger values of X than 
expected for a x 2 with 1 DOF. This can be due to the 
small-sample bias of maximum likelihood methods, 
5 resulting in X if resulting in X x statistics that are not 
X 2 with 1 DOF even if u xi =u yi , for all i. 

We chose A c =25.7, the value at which less than 
0.1% of genes (approximately 6 out of 6200) would be in 
the critical region in the control experiment. This 
10 value was then applied to select differentially-expressed 
genes from the YPR versus YPRG data. 

Figures 2A and 2B show scatter plots of 
estimated ]i y versus \x x values for each gene for the 
control experiment and the YPR versus YPRG experiment, 

15 respectively. The most highly significant genes out of a 
total of 555 selected as significant are shown in 
Table 1. The values shown in Table 1 are in good 
agreement with previous experimental evidence with the 
galactose-induction pathway structural genes GAL1, GAL7 

20 and GAL10 appearing as the top three most significant 
differentially-expressed genes. 



Table 1. Genes Differentially Expressed Between 
Galactose Non-Inducing (YPR) and Inducing (YPRG) 
Conditions . 



10 



GENE 


ROLE 


A 




1 Jf 




GAL1 


galactose metabolism 


95.4 


145 


110644 


766 


GAL 10 


galactose metabolism 


88.1 


109 


36656 


338 


Gal7 


galactose metabolism 


86.7 


59 


76849 


1300 


YNL194C 


unknown 


75.0 


18533 


1360 


0.073 


JEN1 


transport 


72.2 


21124 


889 


0.042 


YNL195C 


unknown 


72.0 


7639 


710 


0.093 


ALD6 


ethanol utilization 


71.5 


9774 


517 


0.053 


RHR2 


glycerol metabolism 


71.1 


1181 


22586 


19 


YMR318C 


unknown 


69.1 


2457 


29930 


12 


HSP26 


diauxic shift 


68.1 


71988 


11435 


0.16 



15 In the scatter plots shown in Figure 2, genes 

with A x >25.7 have significantly different u y and \i x and 
are shown in red. To show detail, axes limits are 
truncated to 45000: the maximum (]i Kf \i y ) observed was 
(1.8xl0 5 , 1.4xl0 5 ) . 

20 Figure 2(C) shows the distribution of four 

(x,y) pairs for two genes in the YPR versus YPRG 
comparison. Samples for each gene are denoted by red or 
black crosses respectively, with corresponding averages 
(<x>, <y>) denoted by squares and MLE-estimated means 

25 (u x ,U y ) denoted by filled circles. Open circles represent 
the estimated means under the added constraint ]U X = y 7 
Pink and gray ellipses define regions containing 95% of 
the error model probability distribution at these 
constrained means for the red and black-colored genes, 
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respectively. Dotted lines of constant ratio, drawn 
through the origin and each constrained and unconstrained 
(U Xf U y ) pair, are shown for reference. In Figure 2C, 
although the genes have similar average expression ratios 
5 <x>/<y> (2.9 versus 3.5 for the red versus black-colored 
gene) , the red-colored gene was significant by the 
likelihood test (A>37.4). The black-colored gene was not 
significant (A=13.8), due to its compatibility with the 
constrained error model. The difference in A arises 
10 because the samples corresponding to the red-colored gene 
are higher in intensity than the samples corresponding to 
the black-colored gene. 

As described in Example I, equation 5 computes 
A for each gene by optimizing the model parameters {\i x 

15 and ji y ) with and without the constraint U x =U y , and 

subsequently compares the likelihood of the (x,y) samples 
under the constrained and unconstrained models. As 
represented by the pink ellipse shown in Figure 2C, the 
four red-colored samples are in the tail of the 

20 probability distribution for the error model with the 
constraint imposed, resulting in a reduced likelihood L 
and thus a relatively high significance value 1. In 
contrast, as represented by the grey ellipse shown in 
Figure 2C, the black-colored samples are relatively well 

25 explained by the constrained error model distribution, 
resulting in a lower value of A. Notably, if the ratio 
statistic r were applied with the commonly-used threshold 
r c =3.0, the black gene would be accepted as significant 
while the red gene would not. 
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Effect of Sample Size on Parameter Estimates 

The more genes and samples per gene are 
available, the more accurate the estimates of |3 and ]a. To 
determine the efficacy of parameter estimation, 
5 representative parameters (3 sim and ]i sim were used to 
randomly simulate data sets of M samples over N genes 
according to the error model equations (2) and (3) 
disclosed in Example I, Values for (3 and y were 
estimated for each data set and the resulting 

10 distribution of |3 over 30 simulations was characterized 
by parameter means <(3> and standard deviations 3 S . In 
simulations with M=50, N=100, parameter estimates were 
tightly distributed around their true values such that 
<(3>=(3 sim ±2% and s 0 <(O.3)<(3> for all parameters g . In 

15 contrast, for very small data sets with M=4, N=100 these 
estimates were highly variable over the 30 simulations 
(s 3 < (0 . 74) <3>) and biased: 3 S im was under- or 
overestimated by 5 percent to 50 percent across the six 
parameters of <(3>. In order to more closely model 

20 experiments performed with a yeast microarray, 
simulations with M=4, N=6000 were also examined. 
Estimates were generally biased but this bias was smaller 
(<3>=3 sim ±25% ) and the variability of estimation also was 
less (s 3 < ( . 05) <(3>) • Thus, with regard to parameter 

25 estimation, a large number of genes appears to at least 
partially compensate for the destabilizing effect of a 
small number of repeats. 

To further study the effect of sample size on 
significance testing in the YPR versus YPRG study, (3, \x 
30 and X were determined using just two of the available 
four samples per gene by drawing one spot per gene over 
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the two replicate hybridizations. In this case, the 
number of genes selected as differentially-expressed was 
less, 227 as compared to 555 using X 1 >25.7, although 85 
percent of these genes were previously identified as 
5 significant when using four samples per gene. The genes 
GAL1, GAL7, and GAL10 also were identified as 
significant, but were no longer among the top ten with 
largest A. While these genes still had a very extreme 
expression ratio (U y /u x ) their intensity samples were by 
10 chance more variable than those of other genes with 

extreme expression ratios and thus their corresponding 
value of X was smaller. 

Ratios of Intensity are Approximately Equal to Ratios of 
Hybridized cDNA 

Although the proposed method identifies genes 
having different mean intensities u xi and iu yi , in order to 
conclude that these genes are differentially-expressed, 
intensity differences or ratios must be at least 
approximately proportional to differences in RNA copy 
number per cell. Since it is expected that either low or 
high copy number could lead to saturation in the measured 
intensity, a series of controlled experiments was 
performed to determine whether this relationship is 
linear over a reasonable range of copy number. 

25 First, a mixture of gal80A cDNA was created by 

extracting mRNA from yeast with a complete deletion of 
the GAL80 gene, which was labelled with Cy3 and Cy5 dyes 
in separate reactions, and subsequently combining the 
reactions into one tube. The mixture was hybridized to a 

30 yeast genome microarray, and the resulting image checked 
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to ensure that intensity was not detectable above 
background for spots representing GAL80 and that all 
spots had roughly equal Cy3 and Cy5 intensities. Next, 
Cy3- and Cy5-labeled DNA sequences corresponding to the 
5 GAL80 open reading frame were added to the gal80A cDNA 
mixture at fixed molar ratios of Cy3:Cy5 dye. 

As shown in Figures 3A and 3B, array 
hybridizations were performed for each of eight 
controlled GAL80 ratios. Data sets consisting of four 

10 (x,y) intensity measurements per gene were obtained at 
each controlled GAL80 ratio by using two spots from a 
forward Cy3:Cy5 labeling scheme and two spots from a 
reverse Cy5:Cy3 labeling scheme. Parameters |3 and ]i were 
determined separately for each data set, and the 

15 corresponding measured ratio for GAL80 was defined as 

iVUx- 

Figure 3C shows a scatter plot of each measured 
ratio versus controlled ratio for the forward-array (red 
dots) or reverse-array (green dots) and demonstrates 

20 that, while saturation occurs at the lower extreme, the 
system is approximately linear over a range of 3 orders 
of magnitude. The ratio of estimated means ia y ,/u K also is 
shown and denoted by open circles. The inset table shows 
the value of X for the GAL80 gene in each of the eight 

25 controlled ratios. The ratio of estimated means U y , /u x is 
denoted by open circles. The inset table shows the value 
of A for the GAL80 gene in each of the eight controlled 
ratios. Except where the controlled ratio was equal to 
one, all measured GAL80 ratios had A> 25.7 and thus were 

30 differentially-expressed by the likelihood test. 
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At the upper end of the investigated range, 
GAL80 was added at 1000 fmol and measured at 32,436 
intensity units as averaged over four samples. Only 14 
genes on the array had higher intensities, the two 
5 largest being TDH3 (81255 units) and EN02 (55766 units) . 
At the lower end of the range, GAL80 was added at 0.2 
fmol and measured at 284 units: approximately 1000 genes 
had lower intensities. These genes are either not 
expressed or are beneath the range of detection. 

The intensities of several genes whose RNA copy 
number per cell has been determined experimentally also 
were determined (Iyer and Struhl, Proc. Natl. Accad. Sci. 
USA 93:5208-5212 (1996)). The RNA corresponding to the 
TRP3 gene has been observed at 1.9 copies per cell in YPR 
media, and had a corresponding average intensity of 597 
(standard deviation of 259) in the YPR condition of the 
YPR versus YPRG array experiment. In contrast, GAL1 mRNA 
is present at less than <0.1 copies per cell in YPR and 
was not significantly above background intensity on our 
yeast array. Thus, most yeast genes, approximately 4000 
to 5000, appear to have intensities within the linear 
range of the microarray system and the lower limit of 
detection is between 0.1 and 1.9 copies/cell. 

Application of the Likelihood Model to Compare and 
25 Contrast Parameters over Different Types of Repeat 
Measurements 

A test microarray having 96 genes spotted 16 
times each was constructed to use the error model to 
compare the combined variability present across an entire 
30 experiment to that introduced during array hybridization 
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and quantitation alone* Ten cultures were grown 
involving identical strains and YPRG conditions, 
independently in separate containers, and RNA prepared 
from each of the ten cultures. Five of the preparations 
5 were labeled using Cy3, while the remaining five were 
labeled using Cy5. The mixtures were combined in Cy3-Cy5 
pairs, and each of the five pairs hybridized to separate 
test arrays. Two types of data sets were drawn from 
these experiments. In the first type of data set, 
10 repeats were drawn from the 16 replicate spots per gene 
on a single array {within- slide data, N=96, M=16) . 

Parameters were estimated by maximum 
likelihood, independently for data sets formed using each 
of the five test arrays. Mean and standard deviation 

15 values over the estimates are shown in Row 1 of Table 2. 
In the second type of data set, repeats were drawn from a 
single spot of each gene on the array over the five 
hybridizations to separate test arrays (between-slide 
data, N=96, M=5) . In this case, parameters 3 were 

20 estimated 16 times, separately for data sets formed using 
each of the 16 spots per gene available on the array (see 
Table 2, Row 2) . Although the multiplicative errors S x 
and £ y have nearly identical standard deviations for the 
within- and between-slide repeats, they are considerably 

25 more correlated within a slide than between slides. In 
addition, the within-slide measurements have less 
variability with regard to the additive error components 
5 K and 5 y . 
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TABLE 2. Comparison of Error Model Parameters for Five 
Within-Slicle and 16 Between-Slide Data Sets. 





Source of 
Variation 






Pe 


o 5x 




5 


within slide 


0.35 


0.306 


0.981 


251 


374 




mean 


(.063) 


(.061) 


(.0069) 


(49) 


(105) 




standard error 












10 


between slides 
mean 

standard error 


0.365 
(.0084) 


0.315 
(.0073) 


0.967 
(.0017) 


422 
(12) 


569 
(13) 



For these optimizations, the parameter p 5 did 
not always converge: it was therefore set to zero during 
parameter estimation and does not appear in Table 2. In 
comparison with other data sets, the prenormalized x r and 
15 y 1 intensities of all 96 genes in the test data were 
moderate to relatively high. Therefore, p 5 was likely 
ill-determined because under the error model, p 6 is 
dominated by p e for larger intensities. 

Although the invention has been described with 
20 reference to the disclosed embodiments, those skilled in 
the art will readily appreciate that the specific 
experiments detailed are only illustrative of the 
invention. It should be understood that various 
modifications can be made without departing from that 
25 spirit of the invention. Accordingly, the invention is 
limited only by the following claims. 



